Precise experimental investigation of eigenmodes in a planar ion crystal 
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The accurate characterization of eigenmodes and eigenfrequencies of two-dimensional ion crystals 
provides the foundation for the use of such structures for quantum simulation purposes. We present a 
combined experimental and theoretical study of two-dimensional ion crystals. We demonstrate that 
standard pseudo-potential theory accurately predicts the positions of the ions and the location of 
structural transitions between different crystal configurations. However, pseudo-potential theory is 
insufficient to determine eigenfrequencies of the two-dimensional ion crystals accurately but shows 
significant deviations from the experimental data obtained from resolved sideband spectroscopy. 
Agreement at the level of 2.5 ■ 10 -3 is found with the full time-dependent Coulomb theory using the 
Floquct-Lyapunov approach and the effect is understood from the dynamics of two-dimensional ion 
crystals in the Paul trap. 

PACS numbers: 03.67.Lx; 45.50. Jf; 37.10.Ty 
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Accurate control of ion crystals is of major importance 
for spectroscopy, quantum simulation or quantum com- 
puting with such experimental platform. Since the in- 
vention of dynamical trapping by Wolfgang Paul [l[ , this 
versatile instrument has been adapted and optimized for 
specific purposes. Charged particles, more specifically 
singly charged ions, are confined in a radio frequency 
(RF) potential, which is formed by tailored electrode 
structures. In the case of the linear Paul trap, one aims 
for a quadrupole field along one z-axis, such that a har- 
monic pseudopotential in x- and y-direction is formed. 
This radial potential strongly confines the ions, while an 
additional weaker axial potential in z-direction is gener- 
ated with static (DC) voltages applied to endcap elec- 
trodes. Trapped ions are cooled by laser radiation Q in 
the potential described by three trap frequencies u x , y ,z 
eventually forming a crystalized structure. 

The conditions of operation are characterized by two 
anisotropy parameters where the radial confinement 
^>(x, y ) typically exceeds the axial DC confinement u) z . 
For sufficiently small values of a(x,y) = ^z/^fxy) ^ nc 
ion crystal is linear and aligned along the weakest, the 
z-trap axis; all ions are placed in the node of the RF 
potential. Spectacular highlights using linear crystals of 
cold ions are the demonstration of quantum logic oper- 
ations d, H[, the generation of entangled states 5, 
sympathetic cooling of ions of different species @, 8], or 
the quantum- logic clock To reach the level of quan- 
tum control, as required in the experiments listed above, 
the first precondition was a complete understanding of 
eigenmodes and eigenfrequencies for such stored linear 

ion crystals [nun. 

For larger numbers of ions, or for larger values of en, 
the linear crystal undergoes a transition to a zigzag struc- 
ture and eventually to a fully crystalline two- or three- 



dimensional structure [H , [l4 1 . Especially interesting are 
planar ion crystals, where usually one of the confining ra- 
dial potentials is much tighter than the axial one, as such 
structures with non-degenerate radial frequencies do not 
rotate Tj| lj| and allow to address and observe individ- 
ual ions. Recent proposals have outlined how to achieve 
laser-induced spin-spin interactions on such 2D ion crys- 
tals and how to use their spatial arrangement for the re- 
alisation of spin lattices exhibiting frustration [171 1 1 81 ] . It 
was also proposed how to study the spin-phonon interac- 
tions coupling the geometric structure of the ion chain to 
the spin-spin interaction of the chain to realize a Peierls 
instability 19| or the Jahn- Teller quantum phase transi- 
tion [lilljTopolo gical defects in the zigzag configura- 



tion were proposed for simulation of quantum effects with 



solitons [22(. It was also proposed, that the dynamics of 



the structural transition from linear to zigzag configura- 
tion may be induced by electr onic ex citations or a fast 
change of the trap parameters 



23-291, and it would al- 



low for the verification of the predicted scaling laws for 
defect formation when traversing the transition, that is 
the Kibble-Zurek mechanism (27], H3- Furthermore, the 
double well structure realized by the two different con- 
figurations of the zigzag ion crystals can be manipulated 
and may thus allow for the creation of a coherent su- 
perposition between these two configurations and hence 
serve as a test bed for decohcrence models 
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Despite these opportunities, no experiments have been 
successful in using planar crystals in a Paul trap. This is 
due to the high complexity of controlling such crystalline 
two-dimensional structures. Required is the knowledge of 
eigenmodes and eigenfrequencies of such crystals, since it 
is important for the cooling and for the design of ion-ion 
interactions which rely on the setting of laser parameters. 
These arc tailored to induce spin-dependent light forces 



2 



via Stark effects 

Here we present a combined experimental and theo- 
retical study of eigenmodes of ions in planar crystals. 
We describe an accurate experimental determination of 
positions of ions and the frequencies of the eigenmodes 
of the crystal structures together with a comparison of 
these data with theoretical expectations based on pseu- 
dopotential approximation (PPT) and a full dynamical 
classical theory for solving the linearized Coulomb prob- 
lem in ion crystals using a Floquet-Lyapunov approach 
(FLT) f~ 



equation is solved 



34 



From observations of multiple sequen- 
tial phase transitions between different structures of the 
ion crystal we map out the phase diagram and find good 
agreement with theoretical predictions following PPT for 
both the transition points and for the positions of indi- 
vidual ions in the different phases. The measurement 
of the eigenfrequencies relies on sideband spectroscopy 
and is applied here to the vibrational frequencies of a 
three-ion zigzag crystal. Surprisingly, some modes ex- 
hibit a significant 37 kHz deviation when compared to 
calculations in PPT. We can understand the measured 
data quantitatively only if the full time dependent solu- 
tion of the trapping potential is taken into account, with 
the eigenfrequencies calculated using FLT. 
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FIG. 1. (color online) Determination of ion positions in planar 
crystals: a) The ion fluorescence near 397 nm is imaged on a 
CCD chip, b) The ion positions (black dots) are determined 
by averaging over 100 exposures and compared with the re- 
sult of a numerical simulation assuming Coulomb repulsion in 
a harmonic trap pseudo-potential (red crosses). The experi- 
mental data allow for a precision of 50 nm as indicated by the 
circle for a la standard deviation. 



The positions of the ions in a crystal are determined 
by the mutual Coulomb repulsion together with a static 
electric potential in the z-direction and a dynamic ra- 
dial trapping potential. A single charged particle in the 
time-dependent potential of a quadrupole Paul trap, or 
the center-of-mass (com) modes of a general crystal of 
ions, obey decoupled, linear Mathieu equations of mo- 
tion in each spatial direction (35[. The RF and DC 
trap voltages, together with the RF frequency fi and the 
mass of the ions m, determine the dimension-free Math- 
ieu parameters in each direction of space i £ {x,y, z} as 
Oj = 4:eUDc/jimfl 2 and g, = 2e[/RF/7,-mf2 2 , where 7, 
and 7j' are geometrical factors denoting the curvature of 
the respective potentials and £ = Qt/2. The Mathieu 



de 



[a -2q . COS (2£)]2, = 0. (1) 

One derives stable solutions with characteristic expo- 
nents Pi w \J a,i + qf /2, and obtains a harmonic time- 
independent pseudo-potential with frequencies Ui = 

Pito/2. 

In the pseudo-potential approximation, the ions ar- 
range in positions corresponding to a stable minimum 
of the time-independent potential, the determination of 
which follows the method for linear crystals (loj , extend- 
ing it to three dimensions. In the time-dependent poten- 
tial, the equivalent of a minimum configuration crystal is 
a periodic solution with the ions oscillating at the RF- 
frequency fi, about well-defined average locations (the 
oscillation denoted as micro- motion). The ion coordi- 
nates in the radial direction obey 



Vn(t) 



- y cos(m) 




(2) 



while in the axial z-direction, the micro-motion is negli- 
gible m. 

In the experiment, the resonance fluorescence of 
Dopplcr cooled 40 Ca + ion crystals near 397 nm is im- 
aged on a CCD camera [3^|. We use a linear Paul trap 
consisting of four cylindrical rods, diameter 2.5 mm and 
at diagonal distance between centres of 4.7 mm as well as 
two DC endcaps at 10 mm distance. For these measure- 
ments we operate with a RF-amplitude of £/rp ~ 300 V pp 
at n/(27r) = 14.62 MHz, and endcap voltages of 350 V, 
yielding secular frequencies u! yz /(2ir)~(316, 111) kHz 
and a much larger uj x . All secular frequencies are deter- 
mined from the resonant electric excitation of the com 
modes, observed as a blurring of the ion crystal im- 
age Q. 

Equilibrium positions are determined from CCD im- 
ages such as shown in Fig. [1] with a 7- ion crystal. Averag- 
ing over 100 exposures and applying a Gaussian fit to the 
data we determine the ion locations, see Fig. [Tp. As the 
planar crystal is observed at an angle of 45°, the y-axis 
is compressed by The magnification of the imaging 
system is determined from the com mode frequency in 
z-direction and the z-distancc in a linear two-ion crystal. 

As the micro-motion amplitude is small compared to 
the diffraction-limited spots of the ion fluorescence, and 
as it consists of a symmetric oscillation about the equi- 
librium positions, the observed fluorescence precisely in- 
dicates the ion positions. The measured data show per- 
fect agreement with the theoretical expectations follow- 
ing PPT at the level of a few parts in 10 5 , see Fig [Tp, 
for a crystal of 7 ions, and similar agreement is found for 
crystals with 6 to 17 ions. 

Structural phase transitions arc induced when the 
value of a is varied. While so far the first lincar-to-zigzag 
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FIG. 2. (color online) Phase transitions for a three-ion (black) 
and a four-ion (red) crystal. Theoretical (PPT) critical a's 
are indicated by vertical dashed lines with the correspond- 
ing crystal configurations. Experimental data are plotted at 
those a where a certain configuration is observed, with an 
error of 0.05 in a. For an ion crystal with even number of 
ions (here N=4) ctsyjn is the relevant parameter, where the 
structure symmetry changes. The measurements are taken in 
a linear micro trap with f2/(27r)=22.7 MHz, Urf ~ 300 V pp 
yielding frequencies W{ y|Z }/(27r) = (0.421, 0.626) MHz and a 
much larger cu x . Inset: Table of calculated critical a. 



transition was investigated [38j , we observe multiple con- 
secutive critical a, see Fig. O The experimentally deter- 
mined positions and the critical values of a agree with 
the prediction of PPT. Comparing the PPT prediction 
and the solution of the full time-dependent equations, 
with q < 0.1, the relative differences for critical a are 
only about 1%, not resolved in the experimental data. 

Normal mode frequencies result from the Coulomb 
forces between ions at their equilibrium positions. For 
the calculation in PPT these forces are expanded in small 
excursions about the equilibrium positions, linearized, 
and the Hessian matrix is solved. Laser spectroscopy pro- 
vides a precise tool to determine eigenfrequencies with a 
relative accuracy of 0.2 percent or better, and we do not 
find agreement between the measurements and the PPT 
calculation. The corresponding values for mode frequen- 
cies for the three-ion crystal under study, in a zigzag con- 
figuration with a w 0.53, are in Tab. 1. Only recently, 
the influence of micro-motion on the frequencies of secu- 
lar modes has been investigated theoretically, which re- 
sults in significant corrections as compared to PPT, even 
for relatively small q- values. 

The breakdown of PPT results from the fast oscilla- 
tions of ions about their equilibrium positions at the RF- 
frequency, as in eq. ([2]), which modifies periodically the 
Coulomb forces between ions. Thus, it is not justified 
to assume static positions for the Hessian matrix. The 
forces between the ions can be expanded in a Fourier se- 
ries, and the resulting equations can be solved in terms 
of decoupled modes [33j , which describe oscillations with 



secular frequencies, superimposed on micro-motion at 
the RF-frcquency. This effect resembles the Lamb shift 
where the Zittcrbewegung of the bound electron leads to 
a modification in the Hydrogen energy levels [3^] . A short 
description of the method of solution, together with code 
for the calculations presented in this Letter, is available 
online 14011. 
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FIG. 3. (color online) Spectrum of vibrational modes in a 
three-ion crystal. The excitation from the S1/2 to D 5 / 2 state 
is plotted as the laser frequency near 729 nm is scanned over 
the resonance between a frequency detuning of ± 2 MHz. We 
observe excitation at the frequencies of red and blue motional 
sidebands, symmetrical around the carrier. The spectrum is 
plotted versus the modulus of frequency detuning, such that 
red and blue sidebands fall on each other. Measurements 
are taken in a linear micro trap with f2/(27r)=35.07 MHz, 
Urf ~ 420 Vpp yielding frequencies w {a . iyiZ) /(27r)={2940(10), 
1695(3), 1238(2)} kHz. We identify zigzag modes near (a) 
1078(2) kHz and (b) 0.714(2) kHz. The blue dashed lines in- 
dicate the expected frequencies for zigzag modes from PPT, 
while the black dashed-dotted lines show the outcome of the 
calculation of FLT. Insets display eigenvectors of zigzag modes 
(a) and (b). The other resonances correspond to mixing fre- 
quencies of two normal mode frequencies. 

We test the FLT prediction experimentally using a 
three-ion crystal and performing resolved sideband spec- 
troscopy on the narrow S\/2 to -D5/2 transition near 
729 nm. In the spectrum (Fig. [3]) wc identify the vibra- 
tional frequencies of a three-ion zigzag crystal. When 
testing PPT, wc use the com modes in axial and ra- 
dial y-direction to generate predictions for the zigzag fre- 
quencies, which deviate from the experimental values by 
37 kHz, and 15 kHz, respectively. For FLT we fit the five 
measurements of the 3 com modes and the 2 lowest pla- 
nar zigzag modes using a weighted-least-squares fit, with 
3 parameters q y , a z and a y (imposing a x = —a z — a y 
as required from the Laplace equation). The weighted- 
least-squares fit norm is a random variable distributed 
like x 2 with 2 degrees of freedom, see . The theoret- 
ical values coincide exactly with both the eigenfrequen- 
cies, fitting the data with about 22% probability with 
all the eigenfrequencies. A similar fitting procedure for 
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Exp. 714(2) 1078(2) 
PPT 730(14) 1041(12) 
FLT 715.1 1078.5 



1238(2) 1695(3) 



1239.5 1690.7 



TABLE I. Frequencies of a three-ion zigzag crystal in units 
of kHz/(2-7r). Experimental uncertainties result from fluctua- 
tions of the trap control voltages, drifts of the laser reference 
cavity or magnetic field noise during the scan. We obtain the 
errors from a statistics of a comparison of red and blue side- 
band frequencies in many scans like that plotted in Fig. 3, 
obtained on the same day and under identical conditions as 
the data. The PPT prediction for the zz-frequencies is calcu- 
lated based on the experimental data for the u) z , y which are 
com modes (whose frequencies are the same in PPT approxi- 
mation and full dynamic theory). Systematic errors including 
the AC-Stark effect < 1 kHz are contained in the error budget. 
From the z- and y-mode frequencies, we first determine the 
ion positions in the pseudo-potential and then the zz-mode 
frequencies. Experimental uncertainties in u z and oj y lead 
to uncertainties for the PPT zigzag mode frequency predic- 
tion of 14 and 12 kHz, respectively. Theory values for FLT 
are obtained as a best fit to all experimentally determined 
frequencies, see text. 



the PPT results in a negligible probability of the order 
10 -11 , reflecting that PPT frequencies do not agree with 
the experimental finding. 
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FIG. 4. (color online) The calculated eigenfrequencies for 
the zigzag modes (a) and (b) in a three-ion crystal plot- 
ted versus the Mathieu parameter q, with a-parameters ad- 
justed such that the com frequencies are constant. The zz- 
eigenfrequencies show a slope with q. The dashed line results 
from the PPT, see Tab. 1, including the prediction error (grey 
shaded) . Experimental data point fit the theoretical FLT ex- 
pectation near q ~ 0.2, where Q./(2n) equals 35.07 MHz. In- 
set: with zoomed part showing the fit of the FLT and the 3cr 
deviation of data against PPT. 

The strong influence of the trap drive on the eigenfre- 
quencies of zigzag modes is further explained by a theo- 
retical analysis, where the Mathieu parameter q is varied 
and the variation of the mode frequencies zz(a) and zz(b) 
is plotted, sec Fig. 2J In the experiment, the q-value is 



fixed and can not varied easily over a large range. The 
data points, near q=0.2, for the zz(a) and zz(b) mode 
hit the FLT prediction and exclude the PPT frequencies. 
The fractional frequency shift reaches for high g-values 
a level of up to 20%. Examining Fig. [5^,, it is under- 
standable why the zz(a) mode is affected the most by the 
micro-motion, and why its frequency increases monoton- 
ically with q. It can be seen from Eq. ([2]) that the micro- 
motion amplitude of each ion is proportional to the nega- 
tive of its y-equilibrium position, so that this motion has 
a large projection on the eigenvector of the zz(b) mode. 
At each point along the periodic trajectory, the restoring 
forces are larger for this mode than at the center, and 
thus its frequency increases with the amplitude of the 
micro-motion, hence with q. 

Conclusion — We calculated and characterized exper- 
imentally the behavior of two-dimensional crystals un- 
der the influence of micro-motion. We show that while 
the pseudo-potential theory fails to explain the exper- 
imental results to an accuracy within 3er-deviation, the 
newly established method predicts the experimental find- 
ing correctly. We could show that FLT provides an ex- 
perimentally verified tool for understanding the static 
and dynamic properties of mesoscopic Coulomb crystals, 
and might prove to be a key ingredient for establishing 
these systems as a new experimental platform. The pre- 
sented results are an essential prerequisite for the success 
of quantum simulation in two-dimensions where the ef- 
fective Hamiltonian depends on the crystal frequencies of 
all normal modes in the presence of micro-motion. Based 
on our results, a theory could be constructed that takes 
into account such effects in the effective Hamiltonian. 

Moreover, the precise knowledge of the behavior of 
the crystal under the influence of micro-motion is cru- 
cial for quantum-to-classical simulation experiments that 
use such transitions in condensed matter or high energy 
effects, such as the Kibble Zurek mechanism. The devia- 
tion of critical a values from the pseudo-potential predic- 
tions is expected to be larger for larger values of q, and 
merits further experimental and theoretical investigation. 
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